Manipulation of organic polyradicals in a single-molecule transistor 
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Inspired by cotunneling spectroscopy of spin-states in a single OPE5-based molecule, we inves- 
tigate the prospects for electric control of magnetism in purely organic molecules contacted in a 
three-terminal geometry. Using the gate electrode, the molecule is reversibly switched between 
three different redox states, with magnetic spectra revealing both ferromagnetic and antiferromag- 
netic exchange couplings on the molecule. These observations are shown to be captured by an 
effective low-energy Heisenberg model, which we substantiate microscopically by a simple valence 
bond description of the molecule. These preliminary findings suggest an interesting route towards 
functionalized all-organic molecular magnetism. 



PACS numbers: 85.65.+h 75.76.+j 85.35.Gv 

I. INTRODUCTION 

Single molecules are of interest as active electronic 
components because of the small size, the possibility 
of self-assembly, and because of their chemical versa- 
tility. As a consequence, the study of single-molecule 
transport!^ is attracting massive interest. Gated three- 
terminal setups (single-molecule transistors are par- 
ticularly useful, since such a geometry allows detailed 
spectroscopy and control of the molecular device 6 . Es- 
pecially interesting is the prospect of exploiting the spin 
degree of freedom in magnetic molecules, either to inter- 
act with and control the charge- or spin currents through 
the system (single-molecule spintronics)PE3] or to store 
and manipulate quantum information (single-molecule 
spin qubits)H In both cases, electric control of the 
spin states is highly desirable since, contrary to mag- 
netic fields, these can be applied locally and varied rela- 
tively fast. Arrays of magnetic atoms, studied with spin- 
polarized STM, were found to be antiferromagnetically 
coupled 1 ^, but ferromagnetic coupling is often more de- 
sirable since a larger net magnetic moment allows easier 
manipulation and readout 14 . 

Here, we discuss the prospect for electrical manipu- 
lation of spin states within a purely organic molecule, 
where the radicals of the charged molecule can be cou- 
pled antiferromagnetically or ferromagnetically depend- 
ing on the relative location of the added charges. This is 
inspired by low-temperature transport measurements on 
a single Coulomb blockaded device, based on a molecular 
"cruciform": thiol end-capped penta phenylene ethyny- 
lene (OPE5) with a vertically disposed redox-active 
tetrathiafulvalene (TTF) unitpl, see Fig. [lja). As we 




FIG. 1: (Color online) (a) Structure of the cruciform 
molecule and a sketch of the measurement setup, (b) Molecu- 
lar spin model of three coupled spins with spin 1/2. Si and S2 
are ferromagnetically coupled, J12 < 0, and both are antifer- 
romagnetically coupled to S3 with J13 « J23 > 0. Increasing 
the gate voltage removes first S3 and then S2. 



show below, the three-terminal transport data reveal two 
ferromagnetically coupled spins, Si and S2, see Fig.[TJb). 
Adjusting the gate voltage to remove an electron from the 
molecule introduces an additional spin, S3, which modi- 
fies the effective interaction of Si and S2 by coupling to 
them antiferromagnetically. This system is analogous to 
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the theoretical suggestion in Refs. I12|16| where Si and 
S2 were used as spin qubits and two-qubit operations 
could be performed using the gate to control their inter- 
action by adding/removing a third spin, S3. Whereas 
the original proposal of Refs. I12|16l involved an intrinsi- 
cally magnetic polyoxometalate molecule, we argue here 
that this can be achieved also in an all-organic molecule 
where a ferromagnetic ground state can be induced by 
gate-controlled oxidation. Thus, the gated junction al- 
lows reversible switching of the molecule between differ- 
ent polyradical states. 



II. METHODS 

The device is made by electromigration 17 of a thin 
gold wire in a solution of the molecules, using a feedback 
mechanism combined with self-breaking to avoid spurious 
gold particle formation in the junctionals]. The current, 
/, was measured as a function of source drain voltage, 
V, and gate voltage, V g , at a temperature T = 2 K. 
The differential conductance, dl/dV, was obtained from 
pointwise linear fits of the /(V) data. 

Due to intentionally low junction fabrication yields to 
avoid formation of junctions comprising several molecules 
contacted in parallel, we only obtained a single device 
with sufficient gate coupling (a = 0.017) and sharp 
enough resonances to allow detailed spectroscopy. In the 
two-dimensional conductance map presented in Fig. [2ja) 
and (b), gate-dependent resonances indicate when V is 
large enough to energetically allow electrons to be added 
to or removed from the molecule by sequential tunnel- 

Below these lines, within the so-called Coulomb dia- 
monds, single-electron transport is suppressed and a fixed 
charge is stabilized by Coulomb blockade. The two cross- 
ing points at V = represent charge degeneracies, and by 
increasing V g going from left to right in Fig. [2|a) and (b) 
a total of two electrons are added to the molecule. Inside 
the Coulomb diamonds transport takes place through co- 
herent many-electron cotunneling processes^, giving rise 
to horizontal conductance features, which can be used 
to e xtrac t information about the molecular spin excita- 
tions 9 10 . In addition, the Kondo effec t! 22 " 24 * gives rise to 
conductance peaks centered at V = 0. 



III. RESULTS AND DISCUSSION 

Figure [2ja) represents a measurement in zero magnetic 
field, while Fig. [2jb) was taken at B = 8T. The dashed 
horizontal lines indicate the excitation energies in each 
charge state. As seen from the line scans in Figs. [2jd)- 
(e), the excitations at finite bias- voltage appear as a split- 
ting of a much broader Kondo-resonance, rather than 
proper inelastic cotunneling steps. Therefore, we choose 
to read off the excitation energies from the ce nter o f the 
peaks rather than from their inflection point s^ 2 ^. For 



the present analysis, we shall merely extract the thresh- 
olds/splittings as an estimate of the exchange couplings 
J12 and J13. 

A. Spin model 

From the clear observation of a Zeeman effect 
(Fig. |2^d)-(h)) we infer that the observed low-energy ex- 
citations must be of magnetic origin. Since changing V g 
across a charge degeneracy point adds one electron at 
the time and we observe no spin blockade 27 , the most 
general effective three-spin model is the one sketched in 
Fig. [TJb) , described by the Hamiltonian 

H = J12S1 • S2 + J13S1 • S3 + J23S2 • S3 

+ <7/X B B.^S;, (1) 

i 

expressed in terms of three distinct 5=1/2 states, mag- 
netic field, B, the Bohr magneton, /ig, and the electron 
^-factor. As we will show below, in the rightmost dia- 
mond S2 = S3 = 0, and in the middle one S3 = 0. 

The model allows an accurate fit of the experimentally 
observed excitation energies and their magnetic field de- 
pendence (compare black full lines and red dashed lines 
in Fig [2ja) and (b)). Importantly, the model also cor- 
rectly predicts which excitations are suppressed due to 
spin-selection rules (see below). The exchange couplings 
and electron ^-factor are fitted to the excitations en- 
ergies in Figs. |2|a) and (b); the best fit is shown by 
the full black lines. We find J 12 = -2.63 ± 0.08 meV, 
J13 ~ J23 = 0.35 ± 0.06 meV and g = 1.37 ± 0.16. We 
ascribe the reduction of the ^-factor to a combination 
of the (Knight-) shift of the spin-resonance frequency due 
to a strong Kondo effectP^ and the aforementioned un- 
certainty in determining the exact excitation energies^. 
The fit is not very sensitive to the exact ratio of the rel- 
atively small couplings J13 and J23, only to their sum. 

In the rightmost Coulomb diamond the strong zero- 
bias Kondo peak and its splitting into peaks at \V\ = 
gfiBB 22 " 2 ^ in magnetic field, indicates a spin- 1/2 ground 
state (Si). The molecular charge must therefore be odd, 
but the absolute charge state cannot be uniquely deter- 
mined from the measurements. We argue below that the 
molecule here is singly reduced, but this is not relevant 
for the present discussion about the spin states. 

Going to the middle diamond, one electron has been 
removed. The presence of a zero-bias Kondo peak still 
implies a spin-degenerate ground state, and removing an 
electron must therefore introduce an additional spin, S2, 
which couples ferromagnetically with Si to form a t riplet 
ground state, resulting in a spin-1 Kondo effect^^. 
That the spin has not changed by more than ±1/2 is con- 
sistent with the absense of spin blockade 27 . Note that, in- 
stead of the usual enhancement of the Kondo peak upon 
approaching the charge degeneracy point a slight depres- 
sion is observed near V g = 2.0 V in Fig. [2] (a). Similar 
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FIG. 2: (Color online) (a) dl / dV measured as a function of V and V g at zero magnetic field, (b) Same as (a), but in a 
magnetic field B — 8T. The red dashed lines represent the excitation energies extracted from the cotunneling spectrum and 
black solid lines are the best fit to the model Hamiltonian |l]). (c) Calculated excitation spectrum [eigenenergies of Eq. |l])] 
as a function of magnetic field in the leftmost diamond, (d) dl/dV as a function of V, obtained by averaging dI/dV(V,V g ) 
between V g = 0.2 V and V g = 0.4 V (leftmost diamond). Black (red) lines are taken at B = (B = 8). (e) and (f) Same as 
(d), but obtained by averaging between V g = 1.4 V and V g = 1.5 V (leftmost diamond) and between V g = 2.4 V and V g = 2.5 V 
(rightmost diamond), respectively, (g) and (h) is magnetic field scan at V g = 1.6 V (triplet ground state) and 2.6 V (doublet 
ground state) respectively. 



observations have been made in Ref. [30], but an expla- 
nation of this requires a more detailed understanding of 
the charge fluctuations in this specific molecule. In a 
magnetic field, the spin-1 Kondo peak splits into low ly- 
ing peaks at \V\ = g^BB. Excitations to the S z = — 1 
component of the triplet are suppressed by spin-selection 
rules and are therefore not present in the experiment. 
The strong inelastic cotunneling lines observed at higher 
bias-voltages indicate transitions to the excited singlet 
configuration and occur at \V\ = — J12 + QI^bB. 

In the leftmost diamond, yet another electron has been 
removed, and the molecule has acquired one more spin, 
(S3). An antiferromagnetic coupling to S12 = Si + S2 
leads to a doublet ground state {S =1/2) and a quadru- 
plet excited state (S = 3/2) as seen in the spectrum 
and magnetic field dependence shown in Fig.^c). This 
explains the merging of the peaks at B = 8T as appar- 
ent from Fig. [5|d); the lowest component of the quadru- 
plet is almost degenerate with the lowest component 
of the doublet, and a new effective doublet comprising 
these two states becomes the groundstate. As a result a 
stronger zero-bias Kondo effect is obtained (compare red 



and black lines at zero bias in Fig. ^d)), meaning; that 
this finite-field doublet is better coupled to the leads than 
the zero-field doublet. A more detailed analysis shows 
that this can indeed be possible, if electrons tunnel via 
both Si and S3. Similar behavior where a Kondo peak 
is split by interaction with another spin, but can be re- 
stored with a magnetic field, has been observed in STM 
studies of coupled magnetic and non-magnetic atoms on 
surfaces^. Unfortunately the device broke down, and as 
a consequence we do not have the full B-field dependence 
in this charge state, but only the conductance map at 
and 8T. The full B-field dependence was however ac- 
quired at the two other charge states (Fig. |2jg-h)). 



The zero-bias Kondo effect arising from the ground 
state doublet at B = (Fig.[2jd), black line) is markedly 
weaker than that observed for the single-spin in the far 
right diamond (Fig. [2^f )) . This can be understood by in- 
spection of the effective three-spin ground state doublet, 
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expressed in terms of individual spin states {SfS^S^) as: 

l ft> = 4 (21 m> - 1 tit) - 1 ;tt» (2) 
I a.) = ±= (i tu) + 1 m> - 21 ut» , (3) 

assuming for simplicity that the two antiferromagnetic 
couplings are equal, i.e. J23 = J13. The amplitudes for 
flipping of this doublet by a simple exchange-cotunnelling 
process at either Si or S3 are therefore reduced by 
Clebsch-Gordan coefficients to respectively 



<ftisrift> 



and I S 3 - 1 ft) 



(4) 



The effective Kondo temperature depends exponentially 
on these amplitudes and is therefore significantly smaller 
than for a single spin. Since the measurement was carried 
out at T < Tk-, a reduced Kondo temperature can easily 
result in a smaller zero-bias conductance. Notice that 
the triplet correlations between Si and S2 makes the ef- 
fective exchange-cotunneling via S3 effectively ferromag- 
netic, thus prohibiting Kondo effect altogether. This sug- 
gests that the strongest tunnel coupling to the molecule 
takes place via Si or S2. 



B. Microscopic model 

Having rationalized the salient features of the data in 
terms of the Hamiltonian in eq. [I] we now proceed to 
give a likely microscopic understanding of this simple 3- 
spin model within a valence bond (VB) picture of the 
molecule^. This gives a description of the low-energy 
excitations of the molecule. The diagram in Fig. [TJa) 
shows one of the possible VB's in the ground state of 
the neutral molecule. This configuration, however is not 
the ground state of the molecule between electrodes. We 
shall use it rather as a neutral, spin singlet, reference 
state. If we first add an electron we arrive at the state 
of the molecule in the rightmost diamond of the stability 
diagram. Adding the electron the following will happen: 
a bond is broken, leaving a (negatively) charged carbon 
atom with no spin, and a neutral carbon with an unpaired 
spin. These two entities (charge and spin) can now move 
through the molecule. The charge will be particularly 
stable at the electrode ends of the molecule due to the 
attraction by image charges^S. Spin and charge will 
not become totally independent, since both disrupt the 
local VB configurations, and staying close minimizes the 
disruption. 

Going to the middle diamond, an electron will be re- 
moved from one of the two TTF-branches of the molecule, 
and the associated unpaired spin will predominantly 
move about in the same branch. Note that the over- 
all molecule is neutral, but we propose that the lowest 
energy configuration has a negative charge close to one 
of the electrodes and a positive charge on one of the sul- 
phur atoms in the TTFs. Once in a while the two spins 




FIG. 3: (Color online) Tri-radical form of the molecular 
cruciform, corresponding to the condition in the leftmost di- 
amond in Fig. [2] One charge (most likely negative), and the 
corresponding spin (Si), is localized close to one end of the 
OPE5 backbone. The vertically disposed TTF unit is oxidized 
twice. The corresponding spins, S2 (S3) are in a meta (ortho)- 
configuration relative to Si, resulting in a ferromagnetic (an- 
tiferromagnetic) coupling. S2 and S3 are para-coupled, re- 
sulting in an antiferromagnetic coupling. 



can perform an exchange process, and this will result in 
either an antiferromagnetic or a ferromagnetic coupling, 
depending on which two branches of the molecule the two 
spins belong to. 

The general rule, which determines the character of the 
magnetic coupling is related to the actual location of the 
involved spins. The carbon atoms in the molecule can 
be divided into two "sublattices" , with nearest neighbor 
atoms belonging to different sublattices. It is a simple 
fact of the VB picture, that to a very good approxima- 
tion, a spin belonging to a given sublattice will stay on 
that sublattice. The rule now is that two spins living 
in the same sublattice will be ferromagnetically coupled, 
while two spins living on different sublattices are antifer- 
romagnetically coupled. 

For the molecule we study, the ground state of the neu- 
tral molecule (the central diamond in the stability dia- 
gram) is a spin triplet. We thus require that the two spins 
live on the same sublattice, which leaves only two possi- 
bilities in the choice of branches with spin and charge. 

Turning to the discussion of the third diamond, we ex- 
tract one more electron. The charge and the spin will 
inhabit a third branch of the molecule. If the first two 
spins were ferromagnetically coupled, the third will nec- 
essarily be antiferromagnetically coupled to the first two 
spins, in accordance with the simplified model used in 
the data analysis. Fig. [3] shows the VB of this state. 

What has been described in words in the previous two 
paragraphs can be put into rigorous mathematics in a 
model, where the VB system can be mapped onto an 
effective Heisenberg model, with as many sites, as there 
are coupled p z orbitals in the molecule^. 
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IV. CONCLUSION 

In conclusion, cotunneling spectroscopy has been used 
to demonstrate the realization of a series of different 
spin states upon successive reduction of a single organic 
molecule. If J13 = J23, the Hamiltonian ([!]) is equivalent 
to the spin-part of Eq. (1) in Ref. \12\ and thus describes 
a two-spin qubit system (Si and S2), with an interac- 
tion which is controllable by the reversible addition / 
removal of S3. However, as long as J13 7^ and J23 ^ 0, 
adding S3, and removing it again after a certain time has 
elapsed, results in a non-trivial change of the two-qubit 
state spanned by Si and S2. Thus, this operation acts 
as a two-qubit gate. 

Using a valence bond description of the molecule we 
have provided a microscopic motivation for the effec- 
tive spin model, which, however, relies on interactions 
with image charges in the electrodes stabilizing a dif- 
ferent charge-configuration than in the molecule in vac- 
uum. Therefore, it cannot be expected that all cruci- 
form molecule devices will behave in the same way as 
the single device investigated here. However, as the 



data stands, lacking the full B dependence in the left- 
most charge state, it demonstrates that the device al- 
lowed electrical control of quantum spins in an organic 
molecule, with the non-trivial spin states being induced 
by the gate- volt age, rather than being an intrinsic prop- 
erty of the neutral molecule. It is an important venue for 
future research to find new molecular structures where 
such controlled redox reactions give rise to the desired 
spin-structure in a more reproducible way. This could 
for example be achieved by adding more intrinsic molec- 
ular redox centers, which can stabilize multiple localized 
charges and spins without the need for interactions with 
the electrodes. 



Acknowledgments 

The research leading to these results has received fund- 
ing from the European Union Seventh Framework Pro- 
gramme (FP7/2007-2013) under agreement no 270369 
(ELFOS) and from FOM. 



1 C. Joachim, J. K. Gimzewski, and A. Aviram, Nature 408, 
541 (2000). 

2 A. Nitzan and M. A. Ratner, Science 300, 1384 (2003). 

3 H. Park, J. Park, A. K. L. Lim, E. H. Anderson, A. P. 
Alivisatos, and P. L. McEuen, Nature 407, 57 (2000). 

4 J. Park, A. N. Pasupathy, J. I. Goldsmith, C. Chang, 
Y. Yaish, J. R. Petta, M. Rinkoski, J. P. Sethna, H. D. 
Abruna, P. L. McEuen, et al., Nature 417, 722 (2002). 

5 S. Kubatkin, A. Danilov, M. Hjort, J. Cornil, J. Bredas, 
N. Stuhr-Hansen, P. Hedegard, and T. Bj0rnholm, Nature 
425, 698 (2003). 

6 E. A. Osorio, K. O'Neill, M. R. Wegewijs, N. Stuhr- 
Hansen, J. Paaske, T. Bj0rnholm, and H. S. J. van der 
Zant, Nano Lett. 7, 3336 (2007). 

7 L. Bogani and W. Wernsdorfer, Nature Materials 7, 179 
(2008). 

8 M. Ruben, A. Landa, E. Lortscher, H. Riel, M. Mayor, 
H. Gorls, H. B. Weber, A. Arnold, and F. Evers, Small 4, 
2229 (2008). 

9 E. A. Osorio, K. Moth-Poulsen, H. S. J. van der Zant, 
J. Paaske, P. Hedegard, K. Flensberg, J. Bendix, and 
T. Bj0rnholm, Nano Lett. 10, 105 (2010). 

10 A. S. Zyazin, J. W. G. van den Berg, E. A. Osorio, 
H. S. J. van der Zant, N. P. Konstantinidis, M. Leijnse, 
M. Wegewijs, F. May, W. Hofstetter, C. Danieli, et al., 
Nano Lett. 10, 3307 (2010). 

11 S. Sanvito, Chem. Soc. Rev. 40, 3336 (2011). 

12 J. Lehmann, A. Gaita-Arino, E. Coronado, and D. Loss, 
Nature Nanotech. 2, 312 (2007). 

13 S. Loth, S. Baumann, C. P. Lutz, D. M. Eigler, and A. J. 
Heinrich, Science 335, 196 (2012). 

14 C. Chappert, A. Fert, and F. N. V. Dau, Nature Mater. 6, 
813 (2007). 

15 K. Jennum, M. Vestergaard, A. Holmen Pedersen, J. Fock, 
J. Jensen, M. Santella, J. J. Led, K. Kilsa, T. Bj0rnholm, 



and M. B. Nielsen, Synthesis 4, 539 (2011). 

16 J. Lehmann, A. Gaita-Arino, E. Coronado, and D. Loss, 
J. Mat. Chem. 19, 1672 (2009). 

17 H. Park, A. K. L. Lim, A. P. Alivisatos, J. Park, and P. L. 
McEuen, Appl. Phys. Lett. 75, 301 (1999). 

18 D. R. Strachan, D. E. Smith, D. E. Johnston, T.-H. Park, 
M. J. Therien, D. A. Bonnell, and A. T. Johnson, Appl. 
Phys. Lett. 86, 043109 (2005). 

19 K. O'Neill, E. A. Osorio, and H. S. J. van der Zant, Appl. 
Phys. Lett. 90, 133109 (2007). 

20 E. A. Osorio, T. Bj0rnholm, J.-M. Lehn, M. Ruben, and 
H. S. J. van der Zant, Jour. Phys. Cond. Mat. 20, 374121 
(2008). 

21 S. DeFranceschi, S. Sasaki, J. M. Elzerman, W. G. van der 
Wiel, S. Tarucha, and L. P. Kouwenhoven, Phys. Rev. Lett. 
86, 878 (2001). 

22 D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, 
D. Abusch-Magder, U. Meirav, and M. A. Kastner, Na- 
ture 391, 156 (1998). 

23 S. M. Cronenwett, T. H. Oosterkamp, and L. P. Kouwen- 
hoven, Science 281, 540 (1998). 

24 J. Nygard, D. H. Cobden, and P. E. Lindelof, Nature 408, 
342 (2000). 

25 A. Weichselbaum, F. Verstraete, U. Schollwock, J. I. Cirac, 
and J. von Delft, Phys. Rev. B 80, 165117 (2009). 

26 S. Schmitt and F. B. Anders, Phys. Rev. Lett. 107, 056801 
(2011). 

27 D. Weinmann, W. Hausler, and B. Kramer, Phys. Rev. 
Lett. 74, 984 (1995). 

28 N. Roch, S. Florens, T. A. Costi, W. Wernsdorfer, and 
F. Balestro, Phys. Rev. Lett. 103, 197202 (2009). 

29 J. J. Parks, A. R. Champagne, T. A. Costi, W. W. Shum, 
A. N. Pasupathy, E. Neuscamman, S. Flores- Torres, P. S. 
Cornaglia, A. A. Aligia, C. A. Balseiro, et al., Science 328, 
1370 (2010). 



6 



L. H. Yu, Z. K. Keane, J. W. Ciszek, L. Cheng, J. M. Tour, 
T. Baruah, M. R. Pederson, and D. Natelson, Phys. Rev. 
Lett. 95, 256803 (2005). 

A. F. Otte, M. Ternes, S. Loth, C. P. Lutz, C. F. Hir- 
jibehedin, and A. J. Heinrich, Phys. Rev. Lett. 103, 107203 
(2009). 



S. Shaik and P. H. Hiberty, A chemist's guide to valence 
bond theory (Wiley, New Jersey, 2008). 
K. Kaasbjerg and K. Flensberg, Nano Lett. 8, 3809 (2008). 
P. Hedegard, to be published. 



